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ABSTRACT 


We investigate the internal and relative motions of the Taurus and Ophiuchus 
star-forming regions using a sample of young stars with accurately measured ra¬ 
dial velocities and proper motions. We hud no evidence for expansion or contrac¬ 
tion of the Taurus complex, but a clear indication for a global rotation, resulting 
in velocity gradients of order 0.1 km s“^ pc“^ across the region. In the case 
of Ophichus more data are needed to reliably establish its internal kinematics. 
Both Taurus and Ophiuchus have a bulk motion relative to the LSR (i.e. a non¬ 
zero mean peculiar velocity) of order 5 km s“^. Interestingly, these velocities are 
roughly equal in magnitude, but nearly exactly opposite in direction. Moving 
back in time, we hud that Taurus and Ophiuchus must have been very near each 
other 20 to 25 Myr ago. This suggests a common origin, possibly related to that 
of Gould’s Belt. 

Subject headings: astrometry — techniques: interferometry — techniques: radial 
velocities — stars: formation — ISM: kinematics and dynamics 
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Introduction 


Ophiuchus (at 120 


2007: 


Torres e 


Wilking et al. 


al. 


2008 




2007, 


Loinard et al. 

2008 

) and Taurus (at 130-160 pc; 

Loinard et al. 


2009 


20 


Kenvon et al 


^ ar e two of the nearest star-forming regions (see 


2008 


, for recent reviews). They have both been 


i nstrnmental in the emergence of our current understanding of low-mass star-formation 


(jShu et ah 


19871) . and have been studied extensively at virt ually all wav e 


instance, there are extensive surveys of Taurus in X-rays by 


by Briceno et a . 


Andrews et al. 


(1993: 1999), near -infrared by 


(120051) . and radio by 


Dzib et al 


observations were o 


(iHaisch et al. 


2002 


20041) . and radio flDzib et al. 


jtained in X-ra.ys (lOzawa et al 


Duchene et al. 


the Spitzer Space Telescope (e.g.. 


Giidel e 


Due 




rene et al 


: al. 


engt 


IS. For 


( 20071) . optical 


(120041) . submillimeter by 


). Similarly, in Ophiuc hus, large-scale 


2005 


20041) . submillimeter 


Gagne et al 


Motte et al. 


1998 


20041). near-infrared 


Johnstone et al 


20131). More recently, both regions have b een targeted by 


Rebull et al. 


20101 


Padgett et al. 


2008) and the Herschel 


Space Observatory (http://www.herschel.fr/cea/gouldbelt/en/). The distance to both 
regions is known very accurately thanks to recent trigonometric parallax me asurements 


20Q7, 

20Q§; 

Torres et al. 

20Q7, 

20Q9, 

2012) 


20121 ). In the case of Taurus, the accuracy of these 


VLBI measurements is sufficient to cha racterize the dept 


reconstruct its 3-dimensional structure ([Torres et al. 


2007 


1 of t 


l e com plex, and crudely 


2009 


20121 ) . 


The proximity of these two regions enables the detection of intrinsically faint sources 
(e.g. substellar objects) and ensures high linear spatial resolution. It also facilitates the 
accurate determination of proper motions, since for a given space velocity, the amplitude 
of the angular displacement diminishes linearly with distance. We will exploit this latter 
property here to derive the three-dimensional velocity vector for a sample of young stellar 
objects distributed across each of the regions. This will be achieved by combining radial 
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velocity measurements from optical and near infrared spectroscopy, with proper motions 
derived from multi-epoch radio interferometric observations -supplemented, of course, by 
the accurate distances mentioned earlier to perform the conversion from angular to space 
velocity. Both conventional interferometers such as NRAO’s Very Large Array (VLA) 
and long baseline interferometers (e.g. NRAO’s Very Long Baseline Array; VLBA) can be 


used to obtain accur ate p roper motions. 


Loinard et al 


found in 


(1200311 and 


Torres et al 


Chandler et al. 


(1200711 . 


sam ples of VLA measurements are shown in 


200511 . while examples of VLBA results can be 


In the present paper, we will collect existing radial velocity and radio proper motion 
measurements available in the literature for young stars in Taurus and Ophiuchus, and 
combine them to construct the 3-dimensional velocity vectors for roughly a dozen young 
stellar systems in Taurus and Ophiuchus. These results will be used to analyze the internal 
and relative kinematics of the Taurus and Ophiuchus star-forming regions. 


2. Compilation of radial velocities and proper motion measurements 


Two important technical points must be made at the outset. The hrst one is that the 


proper motions measured using radio interferometers are, by cons 
reference frame associated with the Solar System barycenter (e.g. 


unction, measured in a 


Thompson et al. 


2nn7h 


The second is relative to the conversion of radial velocities from the LSR to heliocentric 
system. While the results of optical spectroscopy are usually reported in the heliocentric 
system, millimeter spectroscopic observations (that we will use for some of the sources) are 
often reported in the LSR system. For consistency with the proper motion measurements, 
we will express all radial velocities in the heliocentric system. The conversion from LSR to 
heliocentric involves the projection of the Sun motion along the line of sight. As we will see 
below, there is some on-going discussion about the true value of the Solar motion. However, 
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as far as we know, all observatories use the same (fairly old) Solar motion determination, 
corresponding to +20 km s“^ toward B1900 equatorial coordinates (18^; +30°) for the 
conversion from heliocentric to LSR velocities. This is the value that we will use to perform 
the conversion when radial velocities in the literature are expressed in the LSR frame. 


2.1. Taurus 


In the case of the Taurus complex, accurate proper motion measurements are available 
for 7 young stellar systems: Hubble 4, HDE 283571, HP Tan G2, V773 Tan, T Tan, L 1551 
IRS5, and DG Tan. We now briefly present and discuss each source invidually, emphasizing 
the distance, proper motion, and radial velocity measurements. A summary of these 
parameters is provided in Table 1. 


Hubble 4 is a weak line T Tauri star of spectral type K7 located in the dark cloud Lynds 
1495. Its trigonometric parallax (zu = 7.53 ± 0.03 mas, corresponding to 132.8 ± 0.5 pc) 
and its proper motion (/XoCosh = 4.30 ± 0.05 mas yr~^; ju^ = -28.9 ± 0.3 mas yr~^) have 


been m easured using multi-epoch VLBA observations by 


Torres et ah 


fcoOTh . 


Nguven et ah 


(120121 ) End evidence that Hubble 4 may be an SB2 spectroscopic binary, and provide 


multi-epoch measurements of the radial velocities of each of the two stars. The average 
value for the radial velocity of the primary is +18.0 km s“^ with a dispersion of order 1.0 
km s“^ (all expressed in the heliocentric syste m). This is reas o nably consistent with the 


Hartmann et ah 


older measurement of +15.0 ± 1.7 km s ^ by 
will adopt +18.0 ± 2.0 km s“^ for the radial velocity of Hubble 4. 


(1198611 . Gonservatively, we 


HDE 283572 (HIP 20388, V987 Tan) is a G2 star also located in the dark cloud Lynds 
1495 with a VLBA parallax zu = 7.78 ± 0.04 mas (128.5 ±0.6 pc; Torres et ah 2007). Its 
proper motion, also measured with the VLBA, is fZaCosS = 8.88 ± 0.06 mas yr“^; fzs = 
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2 6.6 ifc 0.1 mas yr 


-1 


by 


Nguven et al. 


Torres et al.l 120071) . Its (heliocentric) radial velocity was measured 


(1201211 t o be +14.2 ± 


±1.5 km s ^ reported by 


Walter et al. 


.0 km s This is consistent with the value ±15.0 


198811 


HP Tau G2 (IRAS 04328-2248) belongs, together with the tight binary HP Tan G3, to a 
hierarchical triple system located on the eastern edge of the Taurus complex. Its proper 


motion, fr om VLBA o 
mas yr“^ flTorres et al, 


Dseryations, is fiaCosS = 13.90 ± 0.06 mas yr /i^ = -15.6 ± 0.3 


2009 1. and its VLBA parallax is w = 6.20 ± 0.03 ma s (161.2 ±0.9 


pc; To rres et al. 2009). Its radial velocity is -1 -16.6 ± 1.7 k m s 


-1 


Walter et al 


according to 


Nguyen et al. 


( 19881) . We will adopt here the 


( 2012 1 and ±17.7 ± 1.8 km s ^ according to 
former of these two very consistent values. 

V773 Tau is a well-stud ied quadruple system in the Lynds 1495 cloud, at 132.8 ± 


2.3 pc (ITorres et al. 


20121) . The primary is a tight spectroscopic binary 


spatially resolved in VLBA observations (IBoden et al 


motion has been very well characterized by 


2007: 


Boden et al. 


Torres et al. 


(I2007h and 


(IWe 


ty 


199511 


2012 ). I ts orb ital 


Torres et al. 


(1201 2h by 


combini ng spectroscopic and astrometric data. Two other young stars orbit that central 


binary (iDuchene et al. 


2003 


recently been modeled by 


Boden et al. 


Boden et al. 


2 OI 2 II . The orbit of the nearest of the two has 


(j2012|l . By combining the absolute positions of the 


primary provided b y VLBA observations with the global orbit modeling of the system. 


Torres et al 


(I 2 OI 2 II estimated the proper motion of the system’s barycenter to be /i^ cos 5 
= 8.3 ± 0.5 mas yr“^; = -23.6 ± 0.5 mas yr“^. The (heliocentric) radial velocity of 

the barycenter of the system was estimated by A.J. Boden (private communication) to be 
±16.38 ± 0.52 km s-L 


T Ta u is an extremely well-studied triple system, located at 147.6 ± 0.6 pc (iLoinard et al.l 


(e.g. 


20071. The optically visi 


Duchene et al. 


Die and classical T Tauri star is orbited by an infrared companion 


20021) called T Tau S (the “original” T Tauri star has now been renamed 
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T Tau N). T Tau S is itself a tight binary syste m, where 


been detected both in the infrared and the radio f) Schaefer 


comp 


20iJ; 


ex orbital mo 


Loinard et al. 


ions have 


2003 


20071). 


Since the separation between T Tau S and T Tau N is large (0.7 arcsec, corresponding to 


about 100 AU), the proper motion of T Tau N can reasonably be us ed as a proxy 


proper motion of the entire system. We adopt the value measured by 


Loinard et al. 


or the 


(120031): 


/icCosh = 12.2 ± 0.6 mas yr = -12.7 ± 0.6 mas yr The (heliocentr i c) rad ial 


velocity of T Tau N was measured to be +19.2 ± 0.4 km s ^ by 


Nguven et al. 


i s hig 


ily compatible with the older value of +19.1+1.2 km s ^ reported by 


lll986h 


ll2012h. This 


Hartmann et al. 


L1551 IRS5 is a protostellar binary located in the eponymous dark cloud Lynds 1551 to 
the south-east of the Taurus complex. Given its proximity to T Tau, we will adopt a similar 
distance for Lynds 1551, albeit with an increased uncertainty: 147 ± 5 pc. L1551 IRS5 is 
composed of two protostars separated by about 0.3 arcseconds (about 45 AU), presumably 
in relative orbit. However, given the fairly large separation between the two protostars, 
the orbital motions are small. The absolute and relat ive astrometry of these sources has 


been studied with the VLA by 


Rodriguez et al. 


(120031) . We will adopt the average of the 


two proper motions as the proper motion 
1.6 mas yr“^; /i^ = -21.2 ± 2.5 mas yr“^. 


or the system a s a w hole: fia cos 6 = 13.2 ± 


Fridlund et al. 


(120021) report on high spectral 


resolution observations of the circumbinary disk surrounding the VLA sources, from which 
a (LSR) value of +6.3 ± 1.0 km s“^ can be estimated for the systemic radial velocity. This 
corresponds to +18.3 ± 1.0 km s“^ in the heliocentric frame. 

DG Tau and DG Tau B are located near one another (they are separated by less 


than 1 arcmin), but they do not form a bound system. On the plane o 


locatec 


(l2012a 


about mid-way between L1495 and HP Tau, so we will follow 


the sky, they a re 


Rodriguez et al. 


) in adopting a distance of 150 ± 5 pc (intermediate between 130 pc for L1495 and 



























































160 pc for HP Tau). While DG Tau is a 


younger Class I protostar flWatson et al 


^6 classical T Tauri star, D G Tau B is a somewhat 


2004; 


Taihmau et ah 


of DG Ta u measured wi th the V 
mas yr“^ flRodriguez et ah 


20101). The proper motion 
A is jja cos 6 = 7.5 ± 0.9 mas yr“^; /is = -19.0 ± 0.9 


2012aj). That of DG Tau B, on the o 
3.8 ± 1.9 mas yr“^; fis = -20.6 ± 3.3 mas yr“^ (IRodriguez et ah 


her ha nd, is /Xq cos 6 = 


2012bl ). These are very 


consistent with one another, and we will adopt their weighted mean for the proper motion 
of the DG Tau region: /i^ cos5 = 6.8 ± 0.8 mas yr“^; fis = -19.1 ± 0.9 mas yr“^. 


by 


The fhelioce n 


Nguven et al 


ric) radial velocity of DG Ta u was measurer 


on the other hand, 


( 2012) and 


Zapata et al. 


6.5 km s ^ by 


Bacciotti et al. 


o be +15.4 ± 1.5 km s 


-1 


f 2002f) . For DG Tau B, 


(120151) hnd an LSR systemic velocity of +6.5 ± 1.0 km 


s“^. This corresponds to +16.3 ± 1.0 km s“^ in the heliocentric system. These different 
measurements are highly consistent with each other, and we will adopt 16.1 ± 1.0 km s“^ 
for the (heliocentric) radial velocity of the DG Tau region. 


2.2. Ophiuchus 


Accurate radio proper motions are available for four sources in Ophiuchus: IRAS 16293- 
2422, YLW 15, SI, and DoAr21. The latter two of these sources have measured VLBI 
parallaxes corresponding to a distance of 120 ± 4 pc (Loinard et al. 2008). We will adopt 
this distance for all 4 sources, briefly discussing the specihc case of IRAS 16293-2422 in its 
dedicated section. 


IRAS 16293—2422 is a multiple Glass 0 protostellar system located in the da rk cloud 


ynds 1689N. An estimate of the distance to IRAS 16293-2422 was provided by 


Imai et al. 


( 20071) who used multi-epoch VLBI observations of water masers to obtain a direct 
measurement of the trigonometric parallax. They obtain zu = 5.61q 5 mas, corresponding 
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to d = 178]'^37 pc. However, more recent VLBA water measurements by S. Dzib (private 
communication) are c onsistent with a shor ter distance of order 120 pc, which is also the 


distance estimated by 


Loinard et ah 


(120081) for the Ophiuchus core. It is important to 


with short active phases (e.g.. 


low-mass star 

brming regions are weak, 1 

Clanssen et al. 

199b; 

Desmnrs et al. 

2009) 


20091). As a consequence. 


parallaxes obtained using water masers in low-mass star-forming regions are less reliable 
than those measured from continuum observations of magneticaly active stars (such as SI 
and DoAr21 as repoted by Loinard et ah 2008). Thus, we will adopt 120 pc for the distance 
to IRAS 16293-2422. 


Both the absolute and the relative proper motions of the three protostars in the 


IRAS 16293-2. 


Chandler et ah 


22 sys tem have been measured using multi-epoch VLA observations by 


(120051 1. Two of these protostars (A2 and B) share similar absolute proper 


motions , while the proper motion o 


Loinard (2002) and 


Chandler et al 


the th ird object (Al) is signihcantly different. Following 
(120051) ■ we adopt the mean proper motion of A2 and B 


for the proper motion of the system as a whole, and ascribe the different value measured 
for Al to a significant contribution from its orbital motion. Thus, the proper motion 
adopted for IRAS 162 93-2422 is /!« cos 5 = -16.2 ± 0.9 mas yr“^; fis = ~7.0 ± 1.1 mas yr“^ 


(1 Chandler et al. 


200511 . 


In interferometric millimeter wavelength observations, IRAS 16293-2422 is resolved 


into two cores: one containing 

the B protostar, and the other containing the Al and 

A2 obiects 

(Mundv et al. 

1992' 

. These two condensations have slightly different radial 

velocities 

Prgensen et al. 

2011 

.). Since protostar B is known from the proper motion 


measurements to move little relative to the center of mass of the system, we will adopt the 
radial velocity of component B as a proxy for that of the entire system. This corresponds 
to Rzsr = +2.7 ± 1.9 km s“^, and is equivalent to V), = -7.7 ± 1.9 km s“^ in the heliocentric 
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system. 


YLW 15 (IRAS 16244-2434, IRS 43) is a binary Class I protostar ( Andre et al. 1993 ) 


located in the dark cloud L1682B, near the Ophiuchus core. The pro per motion o : 


members of YLW 15 have been measured from VLA observations by 


Curiel et al. 


both 


20031 


who showed further that source VLAl is the primary of the system, whereas VLA2 is a 
lower-mass companion. Thus, we will adopt the proper motion of VLAl as a proxy for the 
proper motion of the entire system: /i^cosd = -1.4 ± 0.5 mas yr“^; /i^ = -20.8 ± 0.6 mas 


yr 


-1 


In the DCO+ maps of Ophiuchus presented by 


Loren et al. 


(jl990l) . YLW15 is embedded 
in the molecular clump F, whose radial velocity is reported as VJsr = +3.7 ± 0.7 km s“^. 
We will adopt this value of YLW 15 itself, which corresponds to -6.5 ± 0.7 km s“^ in the 
heliocentric reference frame. The assumption that the radial velocity of the molecular gas 
surrounding the star can be taken as a proxy of the radial velocity of the star itself is 
supported by the results of Loinard et al. (2008) and Torres et al. (2009, 2012). 


SI (IRAS 16235-2416, ROX 14, YLW 36) is located in the Ophiuchus core (Lynds 1688). 


It is a B4 star with a mass of about 6 Mq, and it is 

among the brig 

'htest red, near-in 

rared 

far-infrared, X-rav and radio sources in the region ( 

Grasdalen et al 

1973 

Fazio et al. 

1976 

Montmerle et al. 1983: Leons et al. 199ll: Loinard et al. 2008). Its proper motion has been 

measured using multi-epoch VLBA observations by 

Loinard et al. 

2008) 

: Ha cos 6 = -3.88 


± 0.87 mas yr /i^ = -31.55 ± 0.69 mas yr 


We did not hnd any direct (photospheric) radial velocit y measureme n 


literature. However, in the DCO’*' observations reported by 


Loren et al. 


for SI in the 


(119901) . SI is 


located on the edge of clump A. The mean radial velocity of the DCO^ emission for clump 
A is 3.5 km s“^ measured in the LSR. The mean width of the DCO^ lines, on the other 
hand, is 1.0 km s“^, so we adopt visr = 3.5 ± 1.0 km s“^ for the radial velocity of this 
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source. This corresponds to K= ”6.7 ± 1.0 km s ^ in the heliocentric reference frame. 


DoAr 21 (V2246 Oph, Haro 1-6, HBC 637, ROX 8, YLW 26) also belongs to the Ophiuchus 


core. It is a ~ 2.2 Mq star of spe ctral type K 


attributed to a circumstellar disk fjjensen et al 


with an infrared excess around 25 /xm 


variable radio source a nd a bright X-ray source (IMontmerle et al. 


200911. It is associated with a strongly 


1985 


Dzib et ahl 


20131) . The proper motion has been measured by 
fiaCosS = -26.47 ± 0.92 mas yr“^; fis = -28.23 ± 0.73 mas yr~^. 


983; 


Feigelson et al. 


Loinard et al. 


( 120081 ): 


derived from optical spectroscopy was provided by 


Massarotti et al. 


The ra dial velocity 


20051 ): Vr = -4.6 ± 


3.3 km s ^ ( h elioce ntric). This is consistent with the value -6 ± 4 km s ^ reported by 


Jensen et al. 


(120091 ). 


A summary of the proper motion and radial velocity measurements detailed above is 
provided as Table 1. For completeness, we also include the proper motions converted to 
Galactic (£, b) coordinates. From that summary, it is clear that the radial velocities are 
typically accurate to about 1 km s“^. The proper motions, on the other hand, typically 
have a 1-dimensional uncertainty of 1 mas yr“^. At the distance of Ophiuchus and Taurus, 
this also corresponds to about 1 km s“^ errors on the tangential velocity. 


3. Analysis 

3.1. Determination of the 3D velocity vectors 

Since our goal here is to analyze the internal and relative motions of Taurus 
and Ophiuchus, we hrst convert the measured proper motions and radial velocities to 
3-dimensional velocity vectors. We will express these vectors in the rectangular {x, y, z) 
coordinate system commonly used for Galactic studies. The origin of the system is at the 
Sun; the (Ox) axis runs along the Sun - Galactic center direction, the positive direction 
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being toward the Galactic center; {Oy) is in the Galactic plane, orthogonal to (Ox), 
with the positive direction in the direction of Galactic rotation; {Oz) is perpendicular to 
the Galactic plane, oriented toward the Galactic North Pole, thereby making (Oxyz) a 
right-handed coordinate system. From the data in Table 1, both the positions {X, Y, Z) 
and the heliocentric velocities {U, V, W) of each of our stars in the {x, y, z) frame can easily 
be computed. They are listed in Table 2. 


Expressing velocities in the heliocentric system is practical from an observational point 
of view, because the dynamics of the Solar System are so well known that heliocentric 
velocities only contain extremely small systematic uncertainty (i.e. the transformation from 
topocentric to heliocentric velocities is very accurate). For the point of view of Galactic 
Dynamics, however, a Sun-based system is not ideal. In particular, for objects in the Solar 
neighborhood, the Local Standard of Rest (LSR) is preferable. The transformation from 
heliocentric to LSR velocities is effected by subtracting the motion of the Sun relative to 
the LSR from the heliocentric velocities. This is straightforward in principle, but introduces 
signihcant errors in practice, because the Solar motion relative to the LSR is somewhat 
uncertain. As we mentioned earlier, velocities reported in the LSR system normally 
assume a Solar motion of -|-20 km s“^ toward B1900 equatorial coordinates (18^;-|-30°). 
More recent determinations, however, suggest signihcantly different values. Until recently. 


the Hipparcos-based determination of 


Dehnen et ah 


(119981) was widely used. For this 


determination, the components of the Solar motion in the rectangular [x, y, z) coordinate 
system introduced earlier are Uq = 10.00 ± 0.36 km s“^, Vq = 5.25 ± 0.62 km s“^, and Wq 
= 7.17 ± 0.38 km s“^. On the basis of a global analysis of high accuracy trigono metric 


paral 


axes to high-mass star-forming regions distributed across the Galactic plane. 


Reid 


(120091 ) argued in favor of a signihcantly larger value of Vq. T his is supported by a r ecent 


re-analysis of stellar kinematics in the Solar neighborhood by 


Schonrich et al. 


( 2010l ) who 


obtained Uq = 11.1 ± 0.7 km s Vq = 12.2 ± 0.47 km s and Wq = 7.25 ± 0.37 km s 


.,-1 
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Here, we will adopt this latter value to transform the {U, V, W) heliocentric velocities of 
the young stars in Taurus and Ophiuchus into {u,v,w) LSR velocities for those stars. The 
results are given in Table 2. 


3.2. Taurus internal kinematics 


Let us now analyze the 3D velocity vectors in the Taurus complex. In the top row of 
Figure 1, we show their projections onto the {Oxy), (Oxz), and {Oyz) planes. When we 
consider the heliocentric velocities (shown as green arrows in the top row of Figure 1), the 
motions appear highly organized as a result of the dominant reflex motion induced by the 
Sun. When the Solar motion is removed (magenta arrows in the top row of Figure 1), the 
motions appear less clearly organized, although there is still a clear remaining bulk motion, 
particularly in the negative {Ox) direction. This bulk motion is shown as a blue arrow in 
the top row of Figure 1, and will be discussed further below. The smaller value and more 
disorganized aspect of the LSR velocities compared with the heliocentric ones evidently 
reflects the fact that the heliocentric velocities are dominated by (minus) the Solar motion 
itself. While the Sun has a 15-20 km s“^ non-circular (i.e. peculiar) velocity component in 
its orbit around the Galactic center, the Taurus complex is on a much more circular orbit, 
as expected for a region of star-formation. 


To characterize the internal kinematics of the stars in the Taurus complex, we now 
compute the difference {5u, Sv, 6w) between the velocity {u, v, w) of each star in Taurus and 
their mean {u,v,w). For reference, the latter is {u,v,w) = (—5.6, —0.6, —2.0) km s“^, and 
corresponds to {U,V,W) = (—16.7,—12.9,—9.2) km s“^ when expresse d in heliocen tric 
veloci ties. This is very similar to the value obtained independently by iBertout et ah 


(120061) from a larger sample of young stars in Taurus with lower accuracy proper motion 
and distance measurements: {U,V,W) = (—15.4,—11.7,—9.9) km s“^. The projections 
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of {Su,6v,Sw) are shown in the bottom row of Figure 1. They have a fairly random 
appearance, with one-dimensional dispersion = 1.3 km s“^, = 2.1 km s“^, = 3.2 

km s“^. The three-dimensional velocity dispersion is a = = 4.1 km s“^. 

To assess quantitatively the relative importance of random and organized motions 
within Taurus, we proceed as follows. Each star is located at a position relative to the 
center of the complex given by the vector r* and moves relative to the complex at a velocity 
(5v*. To each position vector r*, we associate the unit vector f* = r*/|r*| which simply 
points from the center to each given star in the complex. We will consider two types of 
organized motions: expansion (or contraction) and rotation. The velocity (5v* of each 
star in the complex (measured relative to the complex itself) should be parallel to f* for 
expansion, and anti-parallel for contraction. Thus, the dot product f* . (5v* should be large 
and positive for expansion, and large but negative for contraction. By the same token, the 
cross product f* x (5v* should be small for expansion and contraction. 

Conversely, for large-scale rotation, we expect the cross product f* x to be large 
and the dot product f* . (5v* to be small. For instance, for circular rotation in a disk-like 
structure, (5v* and f* would be orthogonal, so the dot product would be zero and the cross 
product would be maximum. For a 3-dimensional structure such as Taurus, the situation 
would be slightly more complex, but one would certainly expect the cross product to be 
large and the dot product to be small. An alternative way of looking at this issue is that 
the quantity f* x hv* is a proxy for the specihc angular momentum of the complex, which 
is expected to be large for rotation, but small for contraction and expansion. 

We calculated the cross and dot products described above for each star in Taurus, and 
took their mean. Notice that both quantities have dimensions of velocity (this was, indeed, 
the reason for using the unit vector f* rather than r* itself, in the dot and cross products). 
Because the dot product is a measure of expansion, while the cross product is a measure of 
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rotation, we will introduce the following definitions: 

'^exp • ^^*7 


Vrot = f* X (5v*. 

Of course, these quantities are not strictly expansion and rotation velocities, but in view of 
our previous discussion, they can be used as proxies for them. 

For Taurus, we obtain Vexp = -0.15 km s“^ and Vdrc = (—1.55,+2.03, —0.02) km 
s“^. The individual values of the dot and cross products are shown in Table 3, and the 
projections of the individual cross product vectors are shown in the bottom row of Figure 
1. The expansion velocity appears very small compared with the velocity dispersion of ~ 

4 km s“^ measured earlier. This results from the fact that the individual dot products are 
alternatively positive and negative (see Table 3), resulting in a small net mean. Thus, in 
the radial direction, the stellar motions appear to be dominated by a random component 
rather than by an organized expansion or contraction pattern. This is correctly reflected by 
the small absolute value of Vexp- 

The situation for rotation is clearly different. We obtain Vrot of about 2 km s“^ 
comparable with the velocity dispersion. Moreover, the individual cross product vectors 
are clearly not randomly oriented. Instead, their components along the {Ox) axis are 
systematically negative, their components along the {Oy) axis are systematically positive, 
while their component along the {Oz) axis are around zero (see Table 3, Figure 1). This 
suggests that the entire Taurus complex is tumbling with a rotation velocity u in the {Oxy) 
plane. Since Taurus is a few tens of pc across and the rotation velocity is a few km s“^, the 
rotation of Taurus induces velocity gradients of order 0.1 km s“^ pc“^ across the complex. 

The relevance of rotation to the equilibrium of the Taurus complex can be estimated 
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by assuming a homogeneous spherical cloud with Vrot = 2 km s ^ at its edge and comparing 


the gravitational energy w ith the rotationa. 
mass of M ~ 3 X 10^ Mq flUngerechts et al. 


ener 


gy. We assume that Taurus has a total 


2007|). 


19871) and a radius of i? ~ 15 pc flGiidel et ah 


The gravitational energy is given by 


^ 3GM2 

h. ~- 

^grav - 


~ —3 X 10^® ergs, 


while the rotational energy is 


E 


~ 5 X 10^^ ergs 


We then conclude that rotation plays a minor role in the virial equilibrium of Taurus. 

Our analysis of the internal kinematics of Taurus is based on high accuracy radial 
velocity, proper motion, and distance measurements of a limited sample of young stellar 


objects. The comparison (mentioned in S ection 3.2) of 


measured here with the determination by 


Bertout et al. 


he mean bulk motion of Taurus 


(120061) based on a much larger 


sample (but with much less accurate astrometric information) shows that our conclusions 
are trustworthy. It will be very interesting to repeat our analysis with larger samples of 
young st ellar objects when t hey become available. For instance, the Gould’s Belt Distances 


Survey flLoinard et al.l 1201 ill will provide parallaxes and proper motion measurements 


similar t o those used here for tens of young stars in Taurus and other regions. The GAIA 


mission (jde Bruiine 


201211 will provide data with similar accuracy at least for YSOs that 


are not too deeply embedded into their parental dusty cocoons. 
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3.3. The Ophiuchus bulk motion 


The previous analysis could be repeated for Ophiuchus, but the results for internal 
kinematics would be quite uncertain, because there are only four stars with accurate proper 
motions, the stars are highly concentrated (within a few pc of each other), and the necessary 
(but poorly justihed for 2 of the 4 stars) assumption that all are at a common distance. We 


will defer this analysis to a fort 


Gon 


d’s Belt Distances Survey 


icoming paper wher e additional astrometri c results from the 


Loinard et ah 


(1201 ll) will be incorporated fjOrtiz-Leon et ah 


2015h . Here, we will merely use the Ophiuchus results to estimate the mean bulk motion 
of the region. We obtain {u,v,w) = (+4.3,—0.9,+2.4) km s“^, which corresponds to 


{U, V, W) = (—6.8, —13.1, —4.8) km s ^ when expressed in heliocentric velocities. 


3.4. The relative motion between Taurus and Ophiuchus 

Taurus and Ophiuchus are fortuitously located almost symmetrically with respect 
to the Sun: Taurus lies at ~ 145 pc in the direction of the Galactic anti-center, at a 
Galactic latitude ~ -15°. Ophiuchus, on the other hand, lies at about 120 pc in the 
direction of the Galactic center, at a Galactic latitude ~ -15°. In the rectangular system 
that we use throughout this paper, the mean position of the stars that we observed in 
Taurus is (—134.9, +16.8, —42.1) pc, while the mean position of the stars in Ophiuchus is 
(+114.2, —13.7, +34.5) pc. These two positions are almost exactly opposite to one another 
in the rectangular frame where the Sun is at the origin. Remarkably, the mean velocity 
of the stars in Taurus that we calculated earlier {vtau = (—5.6, —0.6, —2.0) km s“^) and 
of those in Ophiuchus {voph = (+4.3, —0.9, +2.4) km s“^) are also almost exactly opposite 
to one another (Figure 2). Both the angle between Vtau and the line joining Ophiuchus to 
Taurus, and the angle between Voph and the line joining Taurus to Ophiuchus are of the 
order of 13° and consistent within the errors with 0°. This strongly suggest a common 
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origin for Taurus and Ophiuchus (Figure 2). Indeed, running the time backwards, we find 
that Taurus and Ophiuchus must have been very near each other about 23.7 Myr ago (this 
is assuming a constant velocity). 

Most young stars and molecular cloud complexes in the Solar neighborhood are 


distributed within an expanding stru cture inc . 


plane, and known as Gould’s Belt (see 


Popp el 


ined by about 15-20° from the Galactic 


19971) for an extensive review). The putative 


center of this structure is located in the G alactic mid-p 


in the direction of the Galactic anti-center (iPerrot et ah 


ane, about 100 pc from the Sun 


20031 ). With this assumed center. 


all nearby substantial star-forming regions in the Solar neighborhood except Taurus can 
be accommodated on an elliptical ring inclined by 17.2° from the Galactic plane, with 
semi- major and minor ax es of 373 and 233 pc, respectively, and a line of node at £ = 


296°. fjPerrot et ahl 120031) . The corresponding dynamical age of the structure is 26.4 Myr, 


remarkably similar to the dynamical age that we derived above for the Taurus-Ophiuchus 
system. Yet, the relation between Taurus and Gould’s Belt is somewhat unclear. Taurus 
appears to be projected in the direction of Gould’s Belt (and is indeed often included in 
Gould’s Belt surveys), at a location intermediate between Perseus and Orion. However, it 
is not contained in the ring tha t dehnes Gould’s B elt. Instead, it is located near the center 


of the Belt (see e.g. Figure 5 in 


Perrot et ah 


20031 


A possible e xplanation for the 


was proposed by 


Olano et ah 


peculiar location of Taur us with respec t to Gould’s Belt 


(119871) ; see also the review by 


Poppel (119971) . In that scheme. 


the Taurus material would have been ejected from a region located somewhere along the 
ring containing the star-forming regions in Gould’s Belt. They argue in favor of a region 
at £ ~ 245° and b ~ -14°. Our analysis of the relative kinematics between Ophiuchus and 


Taurus would be inconsistent with this original position (which would instead 
fairly close to the current position of the Sun), but the mechanism proposed by 


rave to lie 


Olano et ah 
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fll987n could still provide the correct theoretical framework for the observations. In this 


scheme, Taurus and Ophiuchus would originate as the result of an energetic event which 
would have occurred roughly simultaneously with (or only a few Myr after) the creation 
of Gould’s Belt and which would have launched interstellar material on opposite ballistic 
trajectories. It is interesting in this respect to consider the energetics. In combination, 
Taurus and Ophiuchus contain about 5 x 10^ Mq in material, and they are both moving at 
about 5 km s“^. This corresponds to a total kinetic energy of 2 x 10'^® ergs, which is only a 
fraction of the kinetic energy output of a typical core-collapse Supernova explosion. 


4. Conclusions 


In this paper, we have combined radial velocity measurements with high accuracy 
proper motion and parallax determinations for a sample of young stars in Taurus and 
Ophiuchus to characterize both their internal kinematics and their relative motion. We 
hnd no evidence for contraction or expansion of the Taurus complex but fairly conclusive 
indications for global rotation in Taurus. These conclusions will be strengthened once 
additional high quality parallaxes and proper motions become available fo r young stars in 


Taur us and Ophiuchus first as part of the Gould’s Belt Distances Survey flLoinard et ah 


201lh and then from the GAIA mission. 


In addition, we measure the relative velocity of Taurus and Ophiuchus and show that 
they are moving away from each other at a velocity of order 5 km s“^. This points to a 
common origin, some 23.7 Myr ago, possibly related to the phenomena that gave birth to 
Gould’s Belt. 


J.L.R., L.L., S.D., G.N.O, and L.F.R. acknowledge the hnancial support of DGAPA, 
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Table 1. Observational data 


Source 

d (pc) 

fict COS 5 (mas yr 

Us (mas yr 1) 

fi^cosb (mas yr 

^ib (mas yr 1) 

Vr (km s 1) 

IRAS 16293-2422 

120.0 ± 4.0 

-16.2 ± 0.9 

-7.0 ± 1.1 

-15.85+1.06 

7.77+1.1 

-7.7 + 1.9 

YLW 15 

120.0 ± 4.0 

-1.4± 0.5 

-20.8 ± 0.8 

-16.55+0.58 

-12.68+0.58 

-6.5 + 0.7 

SI 

120.0 ± 4.0 

-3.88 ± 0.87 

-31.55 ± 0.69 

-26.25+0.81 

-17.95+0.78 

-6.7 + 1.0 

DoAr 21 

120.0 ± 4.0 

-26.47 ± 0.92 

-28.23 ± 0.73 

-38.69+0.86 

0.85+0.84 

-4.6 + 3.3 

Hubble 4 

132.8 ± 0.5 

-1-4.30 ± 0.05 

-28.9 ± 0.3 

23.95+0.25 

-16.75+0.24 

+ 18.0 + 2.0 

HDE 283572 

128.5 ± 0.6 

-1-8.88 ± 0.06 

-26.6 ± 0.1 

25.53+0.09 

-11.61+0.09 

+ 14.2 + 1.0 

HP Tau G2 

161.2 ± 0.9 

-1-13.90 ± 0.06 

-15.6 ± 0.3 

20.89+0.25 

0.74+0.21 

+ 16.6 + 1.7 

V 773 Tau 

132.8 ± 2.3 

-1-8.3 ± 0.5 

-23.6 ± 0.5 

22.72+0.51 

-10.48+0.53 

+ 16.32 + 0.52 

T Tau N 

146.7 ± 0.6 

-1-12.2 ± 0.6 

-12.7 ± 0.6 

17.59+0.64 

0.95+0.63 

+ 19.2 + 0.4 

L1551 IRS5 

147.0 ± 5.0 

+13.2 ± 1.6 

-21.2 ± 2.5 

24.78+2.18 

-3.12+1.91 

+18.3 + 1.0 

DG Tau (A+B) 

150.0 ± 5.0 

+6.8 ± 0.8 

-19.1 ± 0.9 

18.79+0.91 

-7.62+0.91 

+ 16.1 + 1.0 


Note. — The proper motions in columns 3 and 4 are expressed in equatorial (a, 5) coordinates, while those in columns 5 and 6 are in 
Galactic (£, b) coordinates. 







Table 2. Derived velocities 


Source 

U 

km 

W 

u 

V 

km 

w 

X 

Y 

pc 

Z 

IRAS 16293-2422 

-9.52 

-8.06 

2.15 

1.59 

4.19 

9.4 

114.85 

-12.21 

32.77 

YLW 15 

-5.3 

-8.84 

-8.77 

5.81 

3.41 

-1.52 

114.22 

-14.01 

34.21 

SI 

-5.22 

-14.42 

-11.72 

5.89 

-2.18 

-4.47 

114.04 

-13.81 

34.89 

DoAr 21 

-7.36 

-21.25 

-0.93 

3.75 

-9.01 

6.33 

113.53 

-14.87 

36.09 

Hubble 4 

-17.17 

-11.98 

-14.98 

-6.07 

0.27 

-7.73 

-125.54 

24.77 

-35.54 

HDE 283572 

-14.55 

-13.1 

-10.52 

-3.45 

-0.86 

-3.27 

-122.01 

22.91 

-33.34 

HP Tau G2 

-17.24 

-14.73 

-4.11 

-6.14 

-2.49 

3.15 

-154.43 

11.54 

-45.1 

V 773 Tau 

-16.5 

-11.17 

-10.94 

-5.4 

1.08 

-3.69 

-124.76 

26.03 

-37.37 

T Tau N 

-18.94 

-11.01 

-6.24 

-7.84 

1.24 

1.02 

-136.7 

9.01 

-52.28 

L1551 IRS5 

-16.77 

-16.97 

-8.32 

-5.67 

-4.73 

-1.07 

-138.13 

2.59 

-50.42 

DG Tau (A+B) 

-15.79 

-11.23 

-9.57 

-4.69 

1.02 

-2.32 

-142.89 

20.48 

-40.52 
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Table 3. Dot and cross products for the sources in Taurus 


Source 

r • V 


r X V 


Hubble 4 

0.23 

-1.59 

4.99 

1.63 

HDE 283572 

-2.49 

-3.72 

3.65 

0.87 

HP Tau G2 

0.92 

-0.36 

2.12 

-0.96 

V773 Tau 

-3.19 

-0.32 

2.18 

-1.6 

T Tau N 

0.69 

-1.65 

1.27 

1.06 

L1551 IRS5 

3.01 

-2.82 

0.21 

0.73 

DG Tau (A+B) 

-0.19 

-0.43 

-0.14 

-1.86 
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X (pc) X (pc) Y (pc) 


Fig. 1.— Top: heliocentric (green) and LSR (magenta) velocities for the sources in Taurus 
expressed in the cartesian coordinate system described in the text. The blue arrow shows 
the mean LSR velocity of the Taurus complex. Bottom: The green arrows show Sv, the 
difference between the velocity of each star and the mean velocity of the Taurus complex. 
The blue arrows show the f* x cross product; the black arrow is the mean of these 
cross products. 
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— Taurus-Ophiuchus line 



Fig. 2.— 3D rendering of the relative positions of the Sun, Taurus, and Ophiuchus; the 
arrows show the LSR bulk velocities of Taurus and Ophiuchus. 
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